Classical simulation versus universality in measurement based quantum computation 
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We investigate for which resource states an efficient classical simulation of measurement based 
quantum computation is possible. We show that the Schmidt-rank width, a measure recently intro- 
duced to assess universality of resource states, plays a crucial role in also this context. We relate 
Schmidt-rank width to the optimal description of states in terms of tree tensor networks and show 
that an efficient classical simulation of measurement based quantum computation is possible for 
all states with logarithmically bounded Schmidt-rank width (with respect to the system size). For 
graph states where the Schmidt-rank width scales in this way, we efficiently construct the optimal 
tree tensor network descriptions, and provide several examples. We highlight parallels in the efficient 
description of complex systems in quantum information theory and graph theory. 



I. INTRODUCTION 

The classical description of many-body quantum sys- 
tems, and the classical simulation of their dynamics, 
is generically a hard problem, due to the exponential 
size of the associated Hilbert space 0, 0. Neverthe- 
less, under certain conditions an efficient description of 
states and/or their evolution is possible. This is, for in- 
stance, demonstrated by the density matrix renormaliza- 
tion group method |3j], which allows one to successfully 
calculate ground states of strongly correlated spin sys- 
tems in one spatial dimension using matrix product states 
4J. In this context, the questions 'For which (families 
of) states does an efficient classical description exist?', 
and 'When is an efficient classical simulation of the evo- 
lution of such states under a given dynamics possible?' 
are naturally of central importance. 

Apart from their practical importance, the above ques- 
tions are directly related to more fundamental issues, in 
particular to the power of quantum computation and the 
identification of the essential properties that give quan- 
tum computers their additional power over classical de- 
vices; this relation to quantum computation will be cen- 
tral in this article. In particular, we will study these 
questions from the point of view of the measurement 
based approach to quantum computing, more specifically 
the model of the one-way quantum computer In 
this model, a highly entangled multi-qubit state, the 2D 
cluster state is processed by performing sequences 
of adaptive single-qubit measurements, thereby realiz- 
ing arbitrary quantum computations. The 2D cluster 
state serves as a universal resource for measurement 
based quantum computation (MQC), in the sense that 
any multi-qubit state can be prepared by performing 
sequences of local operations on a sufficiently large 2D 
cluster state. 

When studying the fundamentals of the one-way 
model, two (related) questions naturally arise, which we 
will consider in the following; first, it is asked which re- 
source states, other than the 2D cluster states, form uni- 
versal resources for MQC; second, one may also consider 



the question whether MQC on a given state can be ef- 
ficiently simulated on a classical computer. Naturally, 
these two issues are closely related, as one expects that 
an efficient classical simulation of MQC performed on 
(efficient) universal resource states is impossible. How- 
ever, it is important to stress that classical simulation 
and non-universality are principally different issues. 

The question of which other resource states are also 
universal has been investigated recently in Ref . , where 
the required entanglement resources enabling universal- 
ity were investigated. In particular, it was proven that 
certain entanglement measures, in particular certain en- 
tanglement width measures, must diverge on any univer- 
sal resource, thus providing necessary conditions for uni- 
versality. 

On the other hand, the issue of classical simulation of 
MQC evidently brings us back to the central introductory 
questions posed above. Results regarding the efficient 
simulation of MQC do exist, and it is e.g. known that any 
MQC implemented on a ID cluster state can be simulated 
efficiently @. More generally, the efficient description of 
quantum states in terms of (tree) tensor networks turns 
out to play an important role in this context 0, ll(| . 

In this article we strengthen the connection between 
classical simulation of MQC and non-universality. Our 
starting point will be the no-go results for universal- 
ity obtained in Ref. Q, stating that the entanglement 
monotones entropic entanglement width and Schmidt- 
rank width must diverge on any universal resource; both 
measures are closely related, and we refer to section Hi Al 
for definitions. We then focus on the Schmidt-rank width 
measure, and prove, as our first main result, that MQC 
can be efficiently simulated on every resource state which 
is ruled out by the above no-go result. More generally, 
we prove that MQC can be simulated efficiently on all 
states where the Schmidt-rank width grows at most log- 
arithmically with the system size |ll|. 

Second, along the way of proving the above results, 
we provide a natural interpretation of the Schmidt-rank 
width measure, as we show that this monotone quanti- 
fies what the optimal description of quantum states is in 
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terms of tree tensor networks; this shows that there is in 
fact a large overlap between the present research and the 
work performed in Ref. ^3] regarding the simulation of 
quantum systems using tree tensor networks. 

As our third main result, we show that the Schmidt- 
rank width (and entanglement width) - these are mea- 
sures which are defined in terms of nontrivial optimiza- 
tion problems - can be computed efficiently for all graph 
states. Moreover, for all graph states where the Schmidt 
rank width grows at most logarithmically with the num- 
ber of qubits, we give efficient constructions of the opti- 
mal tree tensor networks describing these states. 

We further remark that the origin of the Schmidt-rank 
width lies in fact in graph theory, and its definition is 
inspired by a graph invariant called rank width. It turns 
out that the study of rank width in graph theory shows 
strong similarities with the study of efficient descriptions 
and simulations of quantum systems, viz. the two intro- 
ductory questions of this article. The similarity is due to 
the fact that, in certain aspects of both quantum infor- 
mation theory and graph theory, one is concerned with 
the efficient description of complex structures in terms 
of tree-like structures. We will comment on the existing 
parallels between these fields. 

Finally, we emphasize that the present work is situ- 
ated in two different dynamic areas of research within 
the field of quantum information theory; the first is the 
study of universality and classical simulation of measure- 
ment based quantum computation, and the second is the 
problem of efficiently describing quantum systems and 
their dynamics. An important aim of this article con- 
sists of bringing together existing results in both fields 
and showing that there is a strong connection between 
them; in particular, we find that the notion of Schmidt- 
rank width has been considered independently in Refs. 
and 01 an d plays an important role in both areas of 
research. In order to establish the connections between 
these two areas in a transparent manner, a substantial 
part of this article is devoted to giving a clear overview 
of which relevant results are known in both fields. 

The paper is organized as follows. In section [H] we 
discuss entanglement width and Schmidt-rank width, 
and their role in universality and classical simulation of 
MQC. In Section UTTl the description of states in terms 
of tree tensor networks is reviewed, and a connection to 
Schmidt-rank width is established. This Section also in- 
cludes our main result, stating that any state with a log- 
arithmically bounded Schmidt-rank width has, in princi- 
ple, an efficient description in terms of a tree tensor net- 
work, and hence any MQC performed on such states can 
be efficiently simulated classically. In Section IIVI these 
results are applied to graph states, and we provide in ad- 
dition an explicit way of obtaining the optimal tree tensor 
network. We discuss the relation between the treatment 
of complex systems in quantum information theory and 
graph theory in section and summarize and conclude 
in section IVll 



II. ENTANGLEMENT WIDTH, 
UNIVERSALITY, AND CLASSICAL 
SIMULATION 

In this section we introduce two related multipar- 
tite entanglement measures called entropic entanglement 
width and Schmidt-rank width and discuss their role in 
the studies of universality of resources for measurement 
based quantum computation (MQC) and in classical sim- 
ulation of MQC. 

These entanglement measures are defined in section 
UTaI In section ITTbI we review the definition of universal 
resources for MQC, and the use of the above measures in 
this study. In section III CI we consider the basic notions 
regarding efficient classical simulation of MQC. Finally, 
in section III Dl we pose the two central questions of this 
article in a precise way; the first question asks about the 
interpretation of the measures entanglement width and 
Schmidt-rank width, and the second deals with the role 
of these measures in the context of classical simulation 
of MQC. 

A. Entanglement width 

The entropic entanglement width -Ewd(IV')) of an 
multi-party state \ip) is an entanglement measure intro- 
duced in Ref. 7]. Qualitatively, this measure computes 
the minimal bipartite entanglement entropy in the state 
\tp), where the minimum is taken over specific classes of 
bipartitions of the system. The precise definition is the 
following. 

Let \ip) be an n-party pure state. A tree is a graph 
with no cycles. Let T be a subcubic tree, which is a tree 
such that every vertex has exactly 1 or 3 incident edges. 
The vertices which are incident with exactly one edge are 
called the leaves of the tree. We consider trees T with 
exactly n leaves V := which are identified 

with the n local Hilbert spaces of the system. Letting 
e = {i, j} be an arbitrary edge of T, we denote by T \ e 
the graph obtained by deleting the edge e from T. The 
graph T\e then consists of exactly two connected compo- 
nents (see Fig. (y, which naturally induce a bipartition 
(A^,B^) of the set V. We denote the bipartite entan- 
glement entropy of with respect to the bipartition 
(A^,,B^) by Ea^,,b^{\^))- The entropic entanglement 
width of the state |?/>) is now defined by 

£wd(|^)) := nun max E A ° T ,B%,(\ip)), (1) 

where the minimization is taken over all subcubic trees 
T with n leaves, which are identified with the n parties 
in the system. 

Thus, for a given tree T we consider the maximum, 
over all edges in T, of the quantity -E J 4f,,Bf,(|V')); tnen 
the minimum, over all subcubic trees T, of such maxima 
is computed. 

Similarly, one may use the Schmidt rank, i.e. the num- 
ber of non-zero Schmidt coefficients, instead of the bipar- 
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FIG. 1: (a) Example of a subcubic tree T with six leaves 
(indicated in blue), (b) Tree T\e obtained by removing edge 
e and induced bipartition (A^,B^). 



tite entropy of entanglement as basic measure. One then 
obtains the Schmidt-rank width, or x~width, denoted by 
XwddV 1 ))- The precise definition is the following. Letting 
Xa^,,b^(\iP}) denote the number of non-zero Schmidt co- 
efficients of \ifj) with respect to a bipartition (Aj.,B!f) 
of V as defined above, the x _ width Q f the state \tp) is 
defined by 

Xwd(|^)) := nrinmaxlogaXA^.B'CI^))- (2) 

It is straightforward to show (cf. Ref. 0) that E w< \ 
is an entanglement monotone [T^. i.e., this measure van- 
ishes on product states, is a local invariant, and decreases 
on average under local operations and classical commu- 
nication (LOCC). The proof can readily be extended to 
Xwd, demonstrating that also x-width is a valid entan- 
glement measure. In fact, using that the Schmidt rank 
is non-increasing under stochastic LOCC, or SLOCC, it 
can be proven that the x _w idth is also non-increasing 
under SLOCC. 

Since the inequality 

log 2 XAs(|V'» >E A , B (m (3) 

holds for any bipartition [A, B) of the system and for any 
state \tp), we have 

Xwd(|<A» > £wd(|^». (4) 

Note, however, that these quantities can show a com- 
pletely different (scaling) behavior. 

It is clear that the definitions of entropic entanglement 
width and Schmidt-rank width are based upon similar 
constructions, where optimizations are performed over 
subcubic trees. Such constructions can of course be re- 
peated for any bipartite entanglement measure; hence 
a whole class of multipartite entanglement measures is 
obtained, which we will call the class of entanglement 
width measures. The entropic entanglement width and 



X-width are two examples of entanglement width mea- 
sures. It would be interesting to consider other examples 
of entanglement width measures, and investigate their 
possible role in quantum information theory. 

The definitions of the above entanglement measures 
are inspired by a graph invariant called rank width, which 
was introduced in Ref. fl3T |. The connection with rank 
width is obtained by evaluating the entropic entangle- 
ment width or x~width in graph states. This is explained 
next. 

First we recall the definition of graph states. Let 
a xi fy, erj denote the Pauli spin matrices. Let G = (V, E) 
be a graph with vertex set V := {1, . . . , n} and edge set 
E. For every vertex a S V, the set N(a) denotes the 
set of neighbors of a, i.e., the collection of all vertices b 
which are connected to a by an edge {a, 6} G E. The 
graph state |G) is then defined to be the unique n-qubit 
state which is the joint eigenstate, with eigenvalues equal 
to 1, of the n commuting correlation operators 

K a :=a^ (g) o-W. (5) 

6eiV(a) 

Standard examples of graph states include the GHZ 
states, and the ID and 2D cluster states, which are ob- 
tained if the underlying graph is a ID chain or a rect- 
angular 2D grid, respectively. We refer to Ref. 0] for 
further details. 

Let r be the n x n adjacency matrix of G, i.e, one has 
= 1 if {a, b} e E and r a & = otherwise. For every 
bipartition (A, B) of the vertex set V, define T(A, B) to 
be the \A\ x \B\ submatrix of T defined by 

T{A,B) :=(T ab ) aeA , beB . (6) 

Using standard graph state techniques it can then be 
shown (see e.g. [l4|) that 

rankp 2 T(A, B) = log 2 x^,s(|G)) 

= E A , B (\G)). (7) 

where rankw 2 X denotes the rank of a matrix X when 
arithmetic is performed over the finite field ¥2 = GF(2). 
Thus, the Schmidt rank and the bipartite entanglement 
entropy w.r.t. any bipartition [A, B) coincide for graph 
states, and are given by the rank of the matrix T(A, B). 
Using the identity J7J, one immediately finds that the 
X-width (and entropic entanglement width) of the graph 
state |G) coincides with the rank width rwd(G) of the 
graph G. The explicit definition of rwd(G) reads [T^ 

rwd(G) := nfinmax ra,nkf 2 T (A^p, B^) (8) 

(where the minimization is again over subcubic trees as 
in the definition of x~width), which, using (JJJ, indeed 
coincides with the x~width of |G). 

Note that the subcubic trees which are considered in 
the definition of rank width are not to be confused with 
the defining graph G of the graph state |G) (the latter can 
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be an arbitrary graph); the subcubic trees merely serve as 
a means of selecting certain bipartitions of the system, in- 
dependent of the state which is considered. For instance, 
if we consider a linear cluster state \Lq) of six qubits, 
corresponding to a graph Lq that is a linear chain, then 
the tree depicted in Fig. ^ corresponds to the optimal 
tree in the definition of the rank- width (and width), 
leading to rwd(L 6 ) = Xwd(\L e }) = 1. 

In section we will further comment on the motiva- 
tions for the definition of rank width, and we will draw 
parallels with the study of complex systems in quantum 
information theory. 

B. Universal resources for MQC 

In Ref. 0| a definition for universality of families of 
states for MQC was put forward, and the use of E^d 
to assess non-universality of states was demonstrated. 
In this section we briefly review the definition and the 
corresponding results. 

Consider an (infinitely large) family of qubit states 

*={|Vi>,|lk>,...}, (9) 

where \ipi) is a state on rrii qubits and rrii < TOj_|_i for 
every i = 1,2,.... This family is called a universal re- 
source for MQC if for each state \<p) on n qubits there 
exists a state \tpi) £ SS? on rrii qubits, with m, > n, such 
that the transformation |^) — > \<p) |0) mi ~™ is possible de- 
terministically by means of LOCC. That is, any state \<j>) 
can be prepared using only states within the family W as 
resource. Equivalently, the action of an arbitrary unitary 
operation U on a product input state |0) n can be imple- 
mented, where now :— U\0) n in the above definition. 
This definition is in the spirit of the model of the one-way 
quantum computer, where sequences of adaptive single 
qubit measurements performed on a sufficiently large 2D 
cluster state allow one to prepare any multi-qubit state. 
The definition of universal resource aims to identify the 
required resources, in terms of entanglement, that allow 
one to perform universal quantum computation in the 
sense specified above. 

In the above definition of universality of a family we 
have not yet considered the efficiency with which states 
can be prepared using members of \P. An efficient uni- 
versal resource ^ is a universal resource having the prop- 
erty that all states that can be efficiently generated with 
a quantum gate network should also be efficiently gener- 
ated from universal resource We refer to Ref. for 
a detailed account on efficient universality. 

In Ref. it was found that any universal resource 
"J - must satisfy the following property. Let E(\(j>)) be 
a functional which is defined on the set of all n-qubit 
states, for all n £ N, and suppose that E(\<j>)) is non- 
increasing under LOCC. More precisely, if \(f>) and \<j>') 
are states on n and n' qubits, respectively, then E(\<j>)) > 
E(\cj)')) whenever the transformation \<j>) — > |0')|O)"~™ is 
possible by means of LOCC. Moreover, let E* denote the 



supremal value of E{\4>)), when the supremum is taken 
over all n-qubit states, for all n £ N (the case E* = oo 
is allowed). Then any universal resource \& must satisfy 
the property 

sup {E(\ip}) | \ip) e *} = E*. (10) 

That is, the supremal value of every entanglement mea- 
sure E must be reached on every universal resource \E'. 
Using the fact that there exist families of quantum states 
where the entropic entanglement width and x~width 
grow unboundedly with the system size (the 2D clus- 
ter states are such examples), it is then straightforward 
to show that any universal family of states ^ must have 
unbounded entropic entanglement width and %— width as 
well. More precisely, one has 0: 

Theorem 1 Let '5 be a universal resource for MQC. 
Then the following statements hold: 

(i) sup {E wd Qip)) | \tp) £ *} = oo; 
(ii) sup {x«.ci(l^» I \*l>) e *} = oo. 

In other words, families VP where the measures E wt ^ or 
Xwd are bounded, cannot be universal. This insight, to- 
gether with the relation between entropic entanglement 
width and %-width and the graph theoretical measure 
rank width, allows one to identify classes of graph states 
as being non-universal since the rank width is bounded 
on such classes. Examples include linear cluster graphs, 
trees, cycle graphs, cographs, graphs locally equivalent to 
trees, graphs of bounded tree-width, graphs of bounded 
clique- width or distance-hereditary graphs. We refer to 
the literature for definitions. 

In the remainder of this paper, we will focus on the 
%-width measure. 

C. Classical simulation of MQC 

Rather than considering the question whether a family 
\& is a universal resource for MQC, one may also con- 
sider the question whether MQC on ^ can be efficiently 
simulated on a classical computer. We will say that ef- 
ficient classical simulation of MQC on a family of states 

is possible, if for every state \ipi) £ $ it is possible to 
simulate every LOCC protocol on a classical computer 
with overhead poly(mj), where rrii denotes the number 
of qubits on which the state \ipi) is defined, as before. We 
remark that an efficient classical description of the initial 
states \tpi) is a necessary, but not necessarily a sufficient 
condition for efficient simulation on a classical computer. 

The issue of classical simulation of MQC has recently 
been considered by several authors. At this point we re- 
mind the reader of what is already known in this context. 
Regarding simulation of MQC on graph states, we recall 
the following results: 

• In Ref. jg it was showed that MQC on ID cluster 
states can be simulated efficiently classically; 
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• In Ref. [T^j ^ was showed that MQC on tree graphs 
can be simulated efficiently classically; 

• In Ref. j2( it was showed that MQC on graphs 
with logarithmically bounded tree width [16| can 
be simulated efficiently classically. 

Note that the above result on tree width implies the two 
other results, as tree graphs (and thus also ID cluster 
graphs) have tree width equal to 1 |l7j . 

More general results, i.e., regarding arbitrary states, 
were obtained in Ref. pjj , where it was shown that MQC 
can be simulated efficiently on all states allowing an ef- 
ficient tree tensor network description. The description 
of quantum systems in terms of tree tensor networks will 
play an important role in the present analysis, and will 
be reviewed in detail in section ITTT1 

Although related, the issues of universality and classi- 
cal simulation in MQC are fundamentally two different 
questions. Most of us expect that any family 4" for which 
classical simulation of MQC is possible, will not be an ef- 
ficient universal resource; this reflects the common belief 
that quantum computers are in some sense exponentially 
more powerful than classical machines - note, however, 
that so far there is no rigorous proof of this statement. 
While one expects the possibility of classical simulation of 
MQC to imply non-universality of a resource 'J, the con- 
verse implication is certainly not believed to hold in gen- 
eral. Indeed, it is highly likely that many non-universal 
families could still be used to implement specific quantum 
algorithms. 

D. Problem formulation 

It is clear that regarding the notion of x _ width, and 
the above issues of universality and classical simulation of 
MQC, a number of open questions remain. In this section 
we formulate two central questions, (Ql) and (Q2), which 
will constitute the main research topics in this article. We 
will first state these questions and then discuss them. 

(Ql) Does there exist a natural interpretation of the \~ 
width measure? 

(Q2) Do there exist resources \& having bounded \~ 
width, which nevertheless do not allow an efficient 
classical simulation of MQC? 

Question (Ql) is concerned with the fact that the defi- 
nition of x _ width seems to be rather arbitrary and not 
intuitive, and solely motivated by the connection to the 
graph theoretical measure rank width. We will, however, 
provide a satisfactory interpretation of this measure in 
the context of quantum information in the next section. 

Question (Q2) is concerned with the question whether 
non-universal resources can still be useful for quan- 
tum computation, in the sense that MQC performed 
on such states is more powerful than classical computa- 
tion. As remarked above, it may well be that there exist 



non-universal families of states where MQC is neverthe- 
less hard to simulate classically. Previous results leave 
open this possibility, as the criteria for non-universality 
and classical simulatability do not coincide. For non- 
universal states detected by the x~width criterion (i.e., 
theorem H (h)) j we will show that this is not the case. 
In section MI CI we will show that MQC can be simu- 
lated efficiently for any family VP which is ruled out by 
the x _ width criterion as being a non-universal resource. 

III. ENTANGLEMENT-WIDTH AND TREE 
TENSOR NETWORKS 

In this section we tackle questions (Ql) and (Q2) as 
stated in the previous section. First we will attach a nat- 
ural interpretation to the %-width measure, as we will 
show that Xwd(IV')) quantifies the complexity of the opti- 
mal tree tensor network (TTN) describing the state l^), 
thus providing a satisfactory answer to question (Ql). 
Moreover, we shall see that this connection with tree 
tensor networks immediately allows us to give a nega- 
tive answer to (Q2): we find that MQC can be simulated 
efficiently on all resources having a bounded %-width. 

These results will be obtained in three main steps. In 
section MI Al we review the notions of tensor networks 
and, more particularly, tree tensor networks. We also 
review results obtained in Ref. , where it was proved 
that LOCC on states specified in terms of efficient TTN 
descriptions can be simulated efficiently; the results in 
Ref. |l(l | will be central ingredients to our analysis. In 
section MI Bl we show how to obtain TTN descriptions 
for arbitrary quantum states. Finally, in section Till Cl wc 
establish the connection between TTNs and x~ width. 

A. Tree tensor networks and efficient simulation of 
quantum systems 

In this section we review the basic notions regarding 
(tree) tensor networks (see also Ref. 0), and the simula- 
tion of quantum systems described by TTNs as obtained 
in Ref. JT^. 

Consider a d\ x • • • x d n complex tensor (l8) 

A := Ai 1 i 2 _.i n , (11) 

where each index i a ranges from 1 to d a , for every a — 
1, . . . ,n. The number of indices n is sometimes called 
the rank of the tensor A. We will call the number D := 
max Q d a the dimension of A. For example, every pure 
7i-qubit state expressed in a local basis, 

l 

\<j>) = J2 A ix-i n \h---in) (12) 
ii,...,i n = 

corresponds to an 2 x • • ■ x 2 tensor of rank ti and dimen- 
sion 2. 
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If A^ and A' 2 ) are two tensors of ranks n\ and ri2, 
respectively, and s and t are integers with 1 < s < n-i 
and 1 < t < ri2, and both the s th index of AW and the 
i th index of A^ 2 ) range from 1 to the same integer d, then 
a sum of the form 
d 

4 (2) mi 

ii---is-i j is+i— in ii— h-i 3 H+i---in y ' 

3=1 

yields a tensor of rank m + n% — 1. This sum is called 
a contraction of the tensors A^ and A' 2 ' . More specif- 
ically, one says that the s th index of A' 1 ' is contracted 
with the t th index of A' 2 ' . A situation where several ten- 
sors are contracted at various indices is 
called a tensor network. The maximal dimension of any 
tensor in the network, is called the dimension of the net- 
work, and will usually be denoted by D in the following. 
Note that every tensor network with n open indices (i.e., 
indices which are not contracted), can be associated in a 
natural way to an n-party pure quantum state. 

We will only consider tensor networks where every in- 
dex appears at most twice in the network. In this case, 
every tensor network can be represented by a graph F in 
the following way. 

• For every tensor a vertex a is drawn. 

• Whenever two tensors A^ and A^are contracted, 
an edge is drawn between the corresponding ver- 
tices a and /3 in the graph. 

• Finally, for every open index of a tensor A^ a \ i.e., 
an index which is not contracted, one draws a new 
vertex and an edge connecting this vertex to the 
vertex a. 

As an example, consider three tensors A^ , A' 2 ' , A^ 
contracted as follows: 

This tensor network has 3 open indices a, 6, c, and the 
indices j, k, I are contracted. The graph underlying this 
tensor network is depicted in Fig. The tensor network 
(|14H is naturally associated with a 3-partite pure state 

w E | E4-US4S I i«>i|b> 2 |c>3, (15) 

abc \ jkl j 

where we introduced local bases {|a)i}, {1^)2}, and {(0)3} 
(the subscripts denote the associated Hilbert spaces of 
the basis vectors). In fact, is an example of a matrix 
product state. Writing 

\^) E A Sli«)i. 

■= E4>> 2 < 
i« - E4 3 ^>s, 



one obtains the shorthand notation 

i^) = E^S ) )i^ ) )i^ ) )- ( 16 ) 
jki 

It is clear that similar shorthand expressions can be ob- 
tained for arbitrary tensor networks. 

A tree tensor network (TTN) is a tensor network where 
the underlying graph is a tree, i.e., a graph with no cycles. 
An example of a TTN is 

V A (1) A {2) A {3) A (4) A {5) A {6) (~\7) 

ijklm 

and the corresponding tree graph is depicted in Fig. [2t>- 
Note that <|14|l is an example of a tensor network which 
is not a TTN. 



(a) (b) 




FIG. 2: Tensor network with three tensors A^ k , A^], A^ 
and three open indices a, b, c corresponding to a cycle graph, 
(b) Tensor network with six tensors A% A\f h , A^ h and 
eight open indices a,b, c, d, e, f, g,h corresponding to a tree 
graph. 

The following definitions regarding TTNs will be im- 
portant below (see theorem^U). Let T be a tree. An open 
edge is an edge which is incident with a leaf of T. An in- 
ner edge is an edge which is not an open edge. Consider a 
TTN with tree T having n open edges, corresponding to 
an n-party state \ip). Let e G T be an inner edge, and let 
(Ay, Bfp) be the corresponding bipartition of the system. 
By partitioning all tensors in the network in two classes 
as induced by the bipartition (Af,, Bfp) and grouping all 
contractions which occur between tensor in the same class 
of the bipartition, one can write the network in the form 

i 

We say that the TTN is in normal form w.r.t the biparti- 
tion (A T , B^) if the vectors {|</>^)} and {l^f,)} are ( U P 
to a normalization) the Schmidt vectors of the state \ip) 
w.r.t the bipartition (A%.,B%,). We say that the TTN is 
in normal form if it is in normal form for all bipartitions 
(Af,, By), where e ranges over all inner edges in T pcj . 

The interest in TTNs in quantum information theory 
lies in the property that the representation of systems 
in terms of TTNs leads to efficient descriptions of states 
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as well as to the possibility of efficiently simulating the 
dynamics of the system. The main results in this context 
were obtained in Refs. @ and [Hj. The latter result will 
be particularly interesting for our purposes, and will be 
reviewed next. 

We will be concerned with TTNs corresponding to sub- 
cubic trees. It can easily be verified that if a TTN cor- 
responds to a subcubic tree, has n open indices, and has 
dimension D, then the TTN depends on at most 0(nD 3 ) 
complex parameters. Therefore, if an n-party state can 
be described by a TTN where D scales at most polynomi- 
ally in n, then can be described by poly(n) complex 
parameters by using this TTN. Hence a family of sys- 
tems allowing an efficient description is obtained. What 
is more, it has been shown that also the processing of 
such systems can efficiently be simulated classically. The 
following result, obtained in Ref. [13, will play an im- 
portant role in the subsequent analysis. 

Theorem 2 If an n-party pure quantum state \ip) 
is specified in terms of a TTN of dimension D, 
where the underlying tree graph is subcubic, then any 
MQC performed on can be classically simulated in 
0(n poly(D)) time. 

Therefore, if D grows at most polynomially with n, then 
the above simulation scheme is efficient. It is noted by 
the authors in Ref. [HJ that there is no restriction in 
considering subcubic trees only, in the sense that any n- 
party state which can be represented by a TTN (with 
arbitrary underlying tree) with poly(n) parameters, can 
also be represented by a subcubic TTN with poly(n) pa- 
rameters. 



B. Description of quantum systems with TTNs 

Theorem shows that, if an efficient TTN description 
is known for a quantum state, then LOCC on this state 
can be simulated efficiently. However, this result does 
not give any information about obtaining an (efficient) 
TTN description of a given state. Note that, if a state 
is specified, there might exist several TTN descriptions, 
some of which might be efficient and some of which might 
not be. In fact, we will see below that, if a subcubic 
tree with n open edges is specified, then any n-party 
state can be represented by a TTN with this specific 
tree structure - although generally tensors of exponential 
dimension in n are required. Therefore, the following two 
questions are naturally raised: 

• If a state and a subcubic tree T are given, what 
is the behavior of the dimension D of the associated 
TTN(s)? 

• If only a state \ip) is given, what is the optimal 
subcubic TTN describing this state, i.e., the one 
with the smallest dimension D? 



Next it is shown that the entanglement in the state 
as measured by the Schmidt-rank, plays a crucial role in 
answering the above questions. We prove the following 
result. 

Theorem 3 Let \tp) be an n-party state and let T be a 
subcubic tree with n leaves which are identified with the n 
parties in the system. Then there exists a TTN descrip- 
tion of with underlying tree T , where the dimension 
D of this TTN is equal to 

log 2 £> = tosix.xa°.,b« (\4>))- (19) 

Moreover, this TTN is in normal form. 

Proof: the proof is constructive. The idea is to stepwise 
compute all tensors associated to the vertices of T, by 
traversing the tree from the leaves to the root, as depicted 
in Fig. |31 First we need some definitions. A vertex of 
T which is not a leaf is called an inner vertex; note that 
every inner vertex has degree 3. We fix one inner vertex 
r and call it the root of the tree T . The depth of a vertex 
is the length of the shortest path from this vertex to 
the root r. We denote by A the maximal depth of any 
inner vertex in T . We refer to Fig. [3] for a schematic 
representation. 



r 




FIG. 3: Subcubic tree with root r, where leaves (correspond- 
ing to the n — 13 parties of the system) are indicated in blue, 
and inner vertices are indicated in red. The tree is arranged 
in such a way that all inner vertices of same depth 6 are on 
the same horizontal line. 

The construction is initialized by considering all inner 
vertices {vi, . . . , vn} of depth A. Every such vertex has 
two open edges, corresponding to two qubits in the sys- 
tem. We let {a a ,b a } be the vertices associated in this 
way to v a , for every a. We then compute all Schmidt 
decompositions w.r.t. the bipartitions ({a a ,b a } - rest of 
the system), i.e., 

M = El*i a) >& (a) >' ( 2 °) 

i 

for every a. The vectors \<j>f^) have support on the qubits 

{a a ,b a }, the vectors \Q ) have support on the rest of 
the system. The Schmidt coefficients are absorbed in the 
latter vectors. 

One then proceeds by computing the tensors associated 
to the inner vertices of depth A — 1, and then to the 
vertices of depth A — 2, ... , up to depth equal to 1, by in 
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every step applying the procedure which will be outlined 
now. 

Let 1 < S < A — 1. For every vertex v, let T v be the 
unique subtree of T such that »eT, and T v is one of the 
two subtrees obtained by deleting the upper edge of v. 
Let T* be the tree obtained by, first, adding one vertex 
v* to T v and connecting v* to v with an edge {v, v*} and, 
second, drawing k open edges at the vertex v*, where k 
is equal to the number of qubits which do not correspond 
to leaves of T v . 

Now, suppose that the following is true: for all inner 
vertices w of depth 5 + 1, a TTN description for \ip) is 
known with tree T^, and all these TTNs are in normal 
form. We then outline a procedure to obtain, for every 
inner vertex v of depth 5, a TTN description for 1-0) with 
tree T*, and all these TTNs are in normal form. 

Procedure. — Consider an inner vertex v of depth 6. 
Let ei, e 2l e 3 denote the edges incident with v, such that 
e\ and e 2 are the lower edges, and e 3 is the upper edge 
as in Fig. 0] Let (Ai, A2, A3) be the unique tripartition 
of the system defined by 

(X U X 2 UX 3 ) := [A%,B%) 
(X 2 ,X 1 UX 3 ) := (A%,B*f) 
(X 3 ,X 1 UX 2 ) := (A%,B$). (21) 

See also Fig. 21 for a simple pictorial definition. 

We then make the distinction between the following 

cases: 

(A) neither ei or e 2 are open edges, i.e., both edges 
connect v to other inner vertices; 

(B) one of these two edges, say e 2 , is an open edge. 

First we consider case (A). Let v\ {v 2 ) be the vertex con- 
nected to v by the edge e\ {e 2 ). By assumption, we have 
TTN descriptions for \tp) with trees T* t and T* 2 which 
are in normal form. Consider these TTN descriptions, 
and group all contractions in such a way that one ob- 
tains Schmidt decompositions of with respect to the 
above bipartitions: 

l^)=El^>l<%>< ( 22 ) 

i=l 

for every a — 1,2, where d a := XXi,Xt (IV 1 )) denote the 
Schmidt ranks, and where Xi denotes the complement of 
Xi (e.g., Xi = X 2 U A3). The Schmidt coefficients have 
been absorbed in the vectors \tp\.)^ Consider also the 
Schmidt decomposition of \tp) w.r.t. the split (A3, Ai U 
X 2 ), using an analogous notation 

d 3 

W=E l^s)^)- (23) 

i=l 

The latter decomposition is not given by TTN so far, 
and has to be calculated separately. Using the above 3 



(a) 




(b) (c) 




FIG. 4: (a) Same subcubic tree as depicted in Fig. |3 
where vertices are re-arranged. We consider an inner vertex 
v of depth 5 = 1 with lower edges ei,e2 and upper edge e^, 
and the corresponding tripartition of the system into groups 
X\, X2X3. (b) Subtree T v and (c) tree (for definition see 
text). 

Schmidt decompositions, we can write 

di 

\i>) = EWOlV^uXa) (24) 
i=l 
dx 

= E Wx^^xM) (25) 

i=l 

= EX^^^kl^uXs) (26) 
«=i j'=i 

= EEEfe)^,)^,^. (27) 
i=i j=i k=i 

where we have used the following arguments and defini- 
tions. In order to go from (12511 to l|26|). we have inserted 
equation 122|) for a = 2 in (|25H : to obtain the last equality 
(|27|l . we have defined the tensor B^ k by 

d 3 

(i>xMkux 3 ) = J2 Bijk \^)- (28) 
fe=i 
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This yields a TTN description of \tp) with underlying 
tree T*. Note that (|27[l implies that the Schmidt vectors 
\^x 1 ux 2 ) are recuperated as 

Wx 1 ux 2 )=J2WxMx 2 ) BlJk > (29) 

i,j 

which shows that the TTN is in normal form w.r.t.the 
bipartition (Xs,Xi U X2). It then immediately follows 
that this TTN is in normal form. This concludes case 
(A). 

Next we consider case (B). Let v\ and V2 be defined 
as above. Note that in this case X2 = {^2}- Con- 
sider again the TTN description and related Schmidt 
decomposition (1221) for a = 1, i.e., for the bipartition 
(Xi,{v2} U A3). Note that the Schmidt decomposition 
for the split ({v 2 }, A1UA3) is not available from the TTN 
since V2 is not an inner vertex, but we will not need it. 
As in (A) , consider also the Schmidt decomposition (|2"51l , 
i.e., for the bipartition (A 3 , X\ U {^2})- We then write 

di 

W) = E WxM^xM) (30) 
1=1 

= EX>Ui^»><^W}> (3 1 ) 

1=1 k=l 

= EE 1^)1^1)1^3) (32) 

i=l k=l 

where we have used the definition 

I^ } >:=<^J^ iU{ m)- ( 33 ) 

This yields a TTN description of with underlying tree 
T* which is again in normal form. This concludes (B). 
This also ends the procedure. 

Note that the assumption of the procedure is trivially 
fulfilled for 6 — A — 1 after the Schmidt decompositions 
(|2L)[I have been computed. The procedure is then applied 
to S = A — 1, A — 2, . . . , 1. After this, all tensors in the 
desired TTN description are known, except the one as- 
sociated to the root r of T. To obtain this final tensor, 
the following steps are taken. Let ei, e2, e 3 be the edges 
incident with r, let v\, i>2, v 3 be the corresponding ver- 
tices of depth 1, and let the tripartition (Xi,X2 1 X 3 ) be 
defined as before. From the previous steps in the algo- 
rithm, we have TTN descriptions for \ip) with trees T* , 
T* 2 and T* 3 which are in normal form. Consider these 
TTN descriptions, and group all contractions as above, 
in such a way that one obtains Schmidt decompositions 
of with respect to the above bipartitions: 

i^> = Ei^>i^>' ( 34 ) 

i=l 



for every a = 1,2,3. A similar derivation as 12411 1271) 
shows that can be written as 

W =i2^i\^x 1 WxMx 3 )B^ k , (35) 
i=i j=i k=i 

where B ljk is defined similarly as above. This expression 
describes as a TTN with tree T, as desired. Moreover, 
it follows from (|35|l that this TTN is in normal form 
w.r.t to the bipartitions (X a ,X a ) for a = 1,2,3. Since 
the TTNs (|34f) were in normal form by construction, this 
implies that the TTN description (|35|l is in normal form 
altogether. 

Finally, it immediately follows that the dimension of 
this TTN is equal to ^ This concludes the proof 
of theorem |21 □ 

Note that theorem|3|proves that, if a subcubic tree with 
n open edges is specified, then any n-party state can be 
represented by a TTN with this specific tree structure. 
The construction presented in the proof of theorem [3] is 
similar to a procedure presented in Ref. of how to 
obtain a matrix product description (which is a particular 
instance of a tensor network) for an arbitrary state | ip); 
there, too, the dimension of the tensor network depends 
on the maximal Schmidt rank of as measured w.r.t a 
specific class of bipartite splits, similar to (but different 
from) cq. ljT§|) . 

C. Connection with \— width 

Theorem |21 now allows us to give a natural interpreta- 
tion of the %-width measure @. Namely, for any state 
one has: 

• XwddV')) is the smallest possible dimension of a 
TTN associated to through the Schmidt de- 
composition construction described in theorem [21 

• the tree T which yields the minimum in @ cor- 
responds to the optimal TTN, i.e., the one with 
smallest dimension 23] . 

These observations fully answer (Ql), the first of the two 
central questions put forward in section III DI of this ar- 
ticle. What is more, we now immediately arrive at a 
satisfactory answer to question (Q2), since theorems 
and |21 (see also Ref. [l(|) imply the following. 

Theorem 4 Let \ip) be an n-party state. Denote \ := 
Xwd{\^P)) } let T be a tree yielding the optimum in the 
definition of \, and suppose that the TTN description of 
with underlying tree T is known. Then any MQC on 
can be simulated classically in 0(n poly(2 x )) time. 

In particular, this result shows that, whenever Xwd is 
bounded on a family of states & = {\ipi),\ip2), . . then 
any MQC on \1/ can be simulated efficiently classically - 
even in linear time in the system size n. This result fully 
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answers question (Q2) in the negative; i.e., the x _ width 
measure, which was originally introduced as a means to 
assess whether a resource \t is universal for MQC, can 
equally well be used to asses whether MQC on \P can 
be efficiently simulated classically. In particular, we have 
found that MQC can be simulated efficiently for any fam- 
ily ^ which is ruled out by the x~width criterion (i.e., 
theorem^] (ii)) as being a non-universal resource. 

Note that theorem 0] even allows one to conclude that 
efficient simulation is possible when Xwd grows at most 
logarithmically with the system size - i.e., it may be un- 
bounded. One observes that if Xwd exhibits this scal- 
ing behavior on a family of states then it is not de- 
tected by the x -w idth universality criterion. This ap- 
parent paradox is resolved by considering the notion of 
efficient universality, which was briefly introduced in sec- 
tion III Bl When this requirement is introduced in the 
definition of universality, the above paradox is resolved 
as follows. One can prove that Xwd (and £" w d) need 
to grow faster than logarithmically with the system size 
on any efficient universal resource. This clearly resolves 
the above apparent contradiction. 

While the above results indeed settle questions (Ql) 
and (Q2), in practical situations one is of course faced 
with the problem whether, when a state \tp) is specified, 
the optimal TTN can be computed efficiently. In particu- 
lar, if theorem^lis to be applied, the following quantities 
need to be computed: 

(a) The quantity x itself; 

(b) an optimal subcubic tree T in the calculation of Xj 

(c) the TTN description of corresponding to the 
tree T. 

It is clear that, for any of the above quantities to be ef- 
ficiently computable, in the least one needs to have an 
efficient description of the state in some form - say, 
a polynomial size quantum circuit leading to the prepa- 
ration of the state, or, in the case where is a graph 
state, the underlying graph or stabilizer description. If 
an efficient description is not available, quantities such 
as e.g. the Schmidt rank w.r.t. some bipartition can 
generally not be computed efficiently, and there is no 
hope of computing e.g. (a) in polynomial time. How- 
ever, it is important to stress that the possibility of an 
efficient description is by no means sufficient to compute 
the quantities (a)-(b)-(c) efficiently. 

Regarding (a) and (b), the optimization in the defi- 
nition of the %-width measure suggests that an explicit 
evaluation of Xwd in a specified state, as well as the deter- 
mination of the optimal subcubic tree, might be a highly 
nontrivial task. However, we note that general results in 
this context are known. In particular, we refer to Ref. 
|13| , where optimization problems of the form 

minmax/(v4y) (36) 

are considered, where / is a function defined on subsets of 
V := {1, ...,n}, / : A C V — > f(A). It has been shown 



that such optimizations can be performed in polynomial 
time in n, i.e., the optimum as well as the tree yielding 
the optimum can be determined efficiently, for a subclass 
of functions / which meet several technical requirements. 
In the next section we will see that the graph states form 
a class of states where these requirements are met, such 
that the calculation of the x-width can be performed ef- 
ficiently. However, the techniques presented in Ref. fl3| 
might be used or generalized to calculate the x-width ef- 
ficiently for classes of states larger than the graph states. 

Regarding (c), it is clear that the optimal TTN de- 
scription of \ip) can only be computed efficiently if this 
TTN description is itself efficient, i.e., if it depends on 
at most poly(n) parameters - this is exactly the case 
when x scales as log(n). If x scales in the latter way, 
then it follows from the procedure outlined in theorem |3 
that the optimal TTN description of \ip) can be obtained 
efficiently given one is able to determine the following 
quantities in poly(n) time: 

(i) the Schmidt coefficients and Schmidt vectors for all 
bipartitions {A^,,B^,), where T is the optimal tree 
in the definition of the x-width. 

(ii) Certain overlaps between Schmidt vectors: in par- 
ticular, the tensor coefficients 

B ijk = (r Xl \^ k x 3 W Xl ux 3 ) (37) 

in eq. (|28() and similar tensors in eq. I|35(l . as well 
as the vectors 

in eq. ffify . 

Thus, a number of conditions need to be fulfilled to ob- 
tain an efficient TTN description, if it exists, for a given 
state. Remarkably, in the next section we show that the 
quantities (a), (b) and (c) can be computed efficiently for 
all graph states. 

As a final remark in this section, note that an efficient 
TTN description (if it exists) of a state \ip) w.r.t. a given 
tree T, can always we obtained efficiently if is already 
specified in terms of an efficient TTN description w.r.t. 
a different tree T 1 . 



IV. GRAPH STATES 

In this section we specialize the results obtained in the 
previous section to graph states. 

A. Simulation of MQC 

Theorem 0] and the connection between x-width of 
graph states and rank width of graphs, allows us to ob- 
tain the following result. 



11 



Theorem 5 Let |G) be a graph state on n qubits. If the 
rank width of G grows at most logarithmically with n, 
then any MQC on \G) can efficiently be simulated classi- 
cally. 

In particular, the above result shows that if rwd(G) is 
bounded on a family Q = {Gi, G 2 , ■ ■ ■ }, then any MQC 
on the set &(Q) = {|Gi), IG2), . . . } can efficiently be sim- 
ulated classically. This provides a complementary result 
to the one obtained in Ref. |3 , where it was proved that 
any family of graphs with bounded rank width cannot 
provide a universal resource for MQC. Therefore, all ex- 
amples given in Ref. of non-universal graph states 
(see also section Hi Bp can also be given here as examples 
of resources on which MQC can be simulated efficiently 
classically. 

Note that theorem [S] supersedes all known results (see 
section 111 C|) on classical simulation of MQC on graph 
states. To see this, let us consider the result in Ref. 
stating that MQC can be simulated efficiently on all 
graph states G with logarithmically bounded tree width 
twd(G). Using the inequality pij ] 

rwd(G) < 4 • twd(G) + 2, (39) 

one finds that, whenever twd(G) scales as log(n) (where n 
is the number of qubits in the system), then also rwd(G) 
scales at most as log(n). Thus, theorem [5] implies that 
MQC can be simulated efficiently on all graph states G 
with logarithmically bounded tree width, and the result 
in Ref. Q is retrieved. This shows that theorem fully 
recovers and generalizes the known results on simulation 
of MQC on graph states. 

Finally, we emphasize that the rank width can be 
bounded on families of graphs which do not at all have 
any tree-like structure, i.e., graphs possibly having many 
cycles; therefore, the presence of cycles in a graph is no 
indication that efficient simulation of MQC on the associ- 
ated state might be hard. One reason of this property is 
that a possible tree structure of a graph does not remain 
invariant under local operations; e.g., the fully connected 
graph and the star graph (one central vertex connected 
to all other vertices) are locally equivalent; the latter is 
a tree graph, the former is not - in fact, the tree width 
of the star graph is equal to 1, whereas the tree width 
of the fully connected graph on n vertices is n — 1 01 • 
Contrary to e.g. the tree-width measure, the rank width 
is a local invariant, thus taking into account such cases. 
Due to these properties, our results prove a significant 
extension to the use of the tree width; indeed, the above 
example unambiguously illustrates the superiority of the 
rank width as a criterion to address the classical simula- 
tion of MQC on graph states. 

B. TTNs for graph states 

In this section we are concerned with the issue whether, 
if a graph state is given, the optimal TTN can be com- 



puted efficiently, i.e., we consider the quantities (a)-(b)- 
(c) as denoted in section Hll CI 

Let G be a graph on n vertices. It was shown in Ref. 
|13| that, for a fixed integer k, the problem "Is the rank 
width of G smaller than k ?" is in the complexity class 
P. Moreover, in Ref. [2^] several polynomial-time so- 
called approximation algorithms for the rank width are 
constructed. When G is given as an input, the (most 
efficient) algorithm either confirms that rwd(G) is larger 
than k, or it outputs a subcubic tree T* such that 

max rank F2 r(A^.,S^,) = 3fc - 1, (40) 

e£T* 

which implies that rwd(G) < 3k — 1. The running time 
of the algorithm is 0(n 3 ). 

These results immediately yield an efficient procedure 
to determine the qualitative behavior of the x~width of 
graph states, and to determine the optimal subcubic tree 
in the calculation of the x~width. More precisely, a possi- 
ble (binary search) approach is the following: first run the 
above algorithm for k — n/2; if the algorithm confirms 
that rwd(G) > n/2, then run the algorithm for k = 3n/4; 
if not, then run the algorithm for k = n/4, etc. This algo- 
rithm is guaranteed to terminate in poly(n) time. After 
the last run of the algorithm, the rank width, and the 
corresponding optimal subcubic tree, is obtained up to a 
factor 3. 

Thus, both quantities (a) and (b) as defined in the dis- 
cussion following theorem^ can be computed efficiently 
for any graph state. 

As for an efficient calculation of quantity (c), we note 
that, for any bipartition of the system, both the Schmidt 
coefficients and the Schmidt vectors can be computed ef- 
ficiently for graph states using the stabilizer formalism; 
moreover, the Schmidt vectors can always be chosen to 
be stabilizer states themselves. This can be proved as fol- 
lows (we only give a sketch of the argument, as it involves 
standard stabilizer techniques). Let |G) be a graph state 
on qubits V := {1, . . . , n}, and let (A, B) be a bipartition 
of V. Let S denote the stabilizer of |G), defined by 

S-=\Y[{K a T« |a! o e{0,l},Vo€vl, (41) 

where the operators K a have been defined in eq. JSJl. 
Thus, S is the commutati ve g roup generated by the op- 
erators K a . One then has [14J 

aGV g£S 

Let Sa be the subgroup of operators in S acting trivially 
on the qubits in V \ A = B. Then 

PA := Tr s |G)(G| = ^ £ g. (43) 
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This operator satisfies 



g&SA H&Sa 



geSA hes A 

The second equality holds since Sa is a group. Denoting 
r := 2I^I|<Sa|~ 1 , it follows that (rp^) 2 = rpA, showing 
that rpA a projection operator. Thus, all nonzero eigen- 
values of this operator are equal to 1 . This shows that all 
nonzero eigenvalues of pA (which are the squares of the 
Schmidt coefficients of \G) w.r.t. the bipartition (A,B)) 
are equal to = 2 - I" 4 I|6>a|. Moreover, as pa has unit 
trace, it follows that 



T&nk(p A ) = 1, 



(45) 



such that the number of nonzero eigenvalues of pa is 
equal to r = ^Sa]' 1 j^. 

The eigenvectors of pa can be computed as follows. 
Let {K A , . . . , Kf} C Sa denote a minimal generating 
set of 5a, where s := log 2 \Sa\- Let {K A +1 , 
additional Pauli operators, chosen in such a way that 



^|l|}be 



K A K A 



(46) 



is a set of commuting and independent operators; such 
a set always exists (though it is non-unique) and can be 
computed efficiently, by using the stabilizer formalism 
(see e.g. ^H). Note that (|46[1 is the generating set of a 
stabilizer state \ifj) on the qubits in A, namely the state 



WW ■■=^A\U( I + K ^- 



(47) 



Moreover, this state is an eigenstate of pa- To see this, 
note that Kf\i/j){ip\ = \i]>){i>\, and thus K A \ip) =\ip), for 
every j = 1, . . . , s. As {K A }j =1 is a generating set of the 
group Sa, this last identity implies that g\ip) — \tp) for 
every g e Sa, and therefore 



(48) 



In order to obtain a basis of eigenvectors, one considers 
the 2l A l- s = 2^l|«S yl |~ 1 stabilizer states \ip as+1 ,. 
with stabilizers generated by 



(49) 



where at — ±1, for every fc = s + 1, . . . , \A\. One can, 
with arguments analogous to above, show that all these 
states are eigenvectors of pa- Moreover, all these states 
are mutually orthogonal; one has 

|^/3„+l,...,/V|) 

= (-l) a * +/J *(^x,...,a |A1 |^. +1 ,...^ l>1 > > , (50) 



for every k 
used that 



s + 1, . . . , where we have respectively 



&S+1, 



(-l) Qfe (^o 



;<X\A\ 



and 



^ fc A IV'/ 33+1 ,...^ |A| > = (-i) /3 '=IV'/ 33+1 ,...,/ 3| , l >. 



(51) 



(52) 



It immediately follows from the identity l|5UI) that the 
states \tpa s+1 ,...,a ]A \ ) are mutually orthogonal. Since there 
are exactly 2^ A \ |<Sa | 1 such vectors, as many as there 
are nonzero Schmidt coefficients, we have computed all 
Schmidt vectors of \G) w.r.t. the bipartition (A, B). Re- 
mark that at this point we only have a stabilizer descrip- 
tion of the Schmidt vectors; if necessary, the expansion of 
these vectors in the computational basis can be computed 
using the results in Ref. |27| . 

This shows that both Schmidt coefficients and Schmidt 
vectors of \G) w.r.t. any bipartition (A, B) can be com- 
puted efficiently, and that the Schmidt vectors can al- 
ways be chosen to be stabilizer states. Moreover, note 
that overlaps between stabilizer states can be computed 
efficiently using stabilizer techniques, and we refer to |28j |. 
where this problem was considered. 

Thus, all necessary ingredients (cf. (i)-(ii) in section 
IIII C(l needed for the efficient construction of the optimal 
TTN of a graph state \G), can be computed efficiently 
when rwd(G) scales as log(n). 

We then arrive at the following result. 

Theorem 6 Let \G) be a graph state on n qubits and 
denote \ '■= Xwd(\G}). Then an optimal subcubic tree T 
in the definition of x can be computed in poly(n) time. 
Moreover, if x scales as log(n) then the TTN description 
of \G) corresponding to T can be computed in poly(n) 
time. 

Note that, in particular, the conditions of the above the- 
orem are fulfilled for all classes of graphs having bounded 
rank width, and thus efficient TTNs can be computed in 
poly(n) time for all such classes. 



C. Example for the cycle graph on n = 6 qubits 

In this section we give an explicit example of the com- 
putation of the rank width, the optimal subcubic tree, 
and the corresponding TTN description of a particular 
graph state, namely the 6-qubit state \C$) associated to 
the cycle graph (or ring graph) Cq on 6 vertices. The 
adjacency matrix T of Cq is the 6x6 matrix 

. 1 ... 1 
1.1... 
■ 1 • 1 • • 
• • 1 • 1 ■ 
... 1 . 1 
1 ... 1 . 



(53) 



where ' • ' denotes an entry equal to zero. 
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1. Rank width and optimal tree 

First we compute the rank width of the graph C§. In 
fact, we will prove that rwd(Cg) = 2. To show this, 
consider the subcubic tree T depicted in Fig. The 
leaves of T are associated to the vertices of C§ in the 
following natural way: first, fix an arbitrary vertex of Cq 
and denote this to be vertex 1; then, starting from vertex 
1, traverse the vertices of Cq in a counterclockwise way, 
and denote the vertices by 2, 3, 4, 5 and 6, respectively; 
these vertices are now associated to the leaves of T by 
identifying vertex 1 with the leftmost leaf of T, vertex 2 
with the second leaf from the left, etc. 

It is now straightforward to show that 

max rank F2 T(A%,,B%) = 2. (54) 

This can be showed by simply computing the ranks of 
all matrices T{A? r ,B^) and picking the largest of these 
ranks. Furthermore, one has 

a T ,(C 6 ) := max rank F2 T(A^,,B^,) > 2 (55) 

eST' 

for every subcubic tree T 1 . This can be seen as follows: 
first, note that <xt'(Cs) > 1 for every T', since 

r&nk V2 T(A,B)>l (56) 

for every bipartition (A,B). Second, suppose that T' is 
a subcubic tree such that aT'(Ce) = 1; we will show that 
this leads to a contradiction. Note that rank F2 r(A, B) 
is equal to 1 if and only if (A, B) is a bipartition of the 
form (one vertex - rest). Moreover, if aT'(Ce) = 1, then 
one must have 

rank F2 T(A C T ,,B^,) = 1 (57) 

for every e G T' . Thus, every bipartition (A^, , Bj,,) must 
be of the form (one vertex - rest); this leads to a contra- 
diction. This shows that the inequality (|55|l is correct. 
We can therefore conclude that 

rwd(C* 6 ) := min a T '(C 6 ) = 2 (58) 

and that the tree T as depicted in Fig. ^ yields the 
optimum. 

At this point we note that here ad hoc methods have 
been used to obtain the above result; however, we remind 
the reader that general algorithms exist to calculate the 
rank width and the optimal tree, as cited in section lTV Al 

2. TTN description 

The computation of the TTN description of |C 6 ) with 
underlying tree T is performed in Appendix 1X1 The re- 
sult is the following: 

abx\X2 ^ abcdx^ cdefx^ efx^XQ ' 

abcdef 

(59) 



where zi,...,Zg S {0,1} and where all indices in the 
sum run from to 1. The pair ab should be regarded as 
one index taking 4 different values, as well as the pairs 
cd and ef. Moreover, one has the following definitions: 

,(1) c r 

Yabx 1 x 2 -~ °a,a;iOb,X2 

/ (2) / \ac+ab+bx 3 +dx 3 

r a bcdx 3 — \ L ) 

^e/lsx« := *e,*B <*/,<■* • ( 60 ) 

V. COMPLEX SYSTEMS VERSUS TREE 
STRUCTURES IN QIT AND GRAPH THEORY 

We have seen that the x _ width of a graph state is equal 
to the rank width of the underlying graph. There is in 
fact a striking parallel between the motivations for the 
definitions of rank width of graphs and of x - width f 
general quantum states, on which we comment here. 

As explained above, the x _w i ( ith gives information 
about the optimal TTN which describes a given quantum 
state. The interest in such TTNs naturally arises due to 
the fact that the dynamics of quantum systems which 
allow TTN descriptions with sufficiently small dimension 
D, can be simulated efficiently on a classical computer. 
These and similar techniques (cf. e.g. the matrix product 
states formalism) are invoked because the efficient classi- 
cal simulation of general quantum systems can be a very 
difficult problem. Thus, in spite of the general hardness 
of this simulation problem, it becomes tractable when re- 
stricted to the class of those systems with efficient TTN 
descriptions. 

In graph theory an analogous situation occurs. While 
many interesting problems are hard to compute on gen- 
eral graphs, they become tractable for those classes of 
graphs which can be associated, through certain con- 
structions, with tree structures. The simplest examples 
are of course the tree graphs themselves, which are in 
some sense the simplest instances of graphs; and indeed, 
many difficult problems become efficiently solvable, or 
even trivial, on trees. However, this is far from the whole 
story. In graph theory one has considered a variety of so- 
called width parameters, which all measure, in different 
ways, how similar a graph is to a tree graph. Examples 
are rank width, tree width, clique-width, path-width, 
and branch-width. It has been shown that for fami- 
lies of graphs where a given width parameter is bounded, 
large classes of (NP-)hard problems have efficient solu- 
tions. For example, the problem of deciding whether a 
graph is 3-colorable, which is a NP-hard, is efficiently 
solvable when restricted to classes of graphs of bounded 
rank width. The graph theoretical results in this context 
are often very general and far-reaching; e.g., it has been 
show that all graph problems which can be formulated in 
terms of a certain mathematical logic calculus, have effi- 
cient solutions when restricted to graphs of bounded rank 
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width. We refer to Ref. [29j for an accessible treatment 
of these and related issues. 

Thus, in certain aspects of both quantum information 
theory and graph theory there is a natural interest in 
using tree structures for the approximation of complex 
systems. Moreover, there seems to be a strong parallel in 
the explicit constructions which are used in both fields. 
A striking example is obtained here, as the rank width 
of graphs exactly coincides with the %-width measure on 
raph states. As a second example, it was found in Ref. 
that the efficient contraction of large tensor network 
is directly related to the tree width of the underlying 
graphs. The present authors believe that the aforemen- 
tioned parallel can significantly be exploited further. 



VI. CONCLUSION 

In this paper we have considered the possibility to clas- 
sically simulate measurement based quantum computa- 
tion. We have shown that all states with a bounded or 
logarithmically growing Schmidt-rank width can in fact 
be described efficiently, and moreover any one-way quan- 
tum computation performed on such states can also be 
simulated efficiently. We have given an interpretation of 
the Schmidt-rank width, a measure that has its origin in 
graph theory, in terms of the optimal tree tensor network 
describing a state. We have also provided a constructive 
procedure how to obtain the optimal TTN, and discussed 
the requirements that this can be done efficiently. For 
graph states, we have explicitly constructed the corre- 
sponding TTN, and provided an efficient algorithm to do 
this for any graph state where the underlying graph has 
bounded or logarithmically growing rank width. These 
results on efficient simulation complement recent find- 
ings on universality of states, in the sense that all states 
that are found to be non-universal resources for MQC 
using the Schmidt-rank width criteria (i.e. which have 
bounded Schmidt-rank width) can also be simulated effi- 
ciently on a classical computer. The connection to com- 
plexity issues in graph theory, also highlighted in this 
paper, seems to provide future possibilities for a fruitful 
interchange of concepts and methods between the fields 
of quantum information and graph theory. 
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APPENDIX A: OPTIMAL TTN DESCRIPTION 

OF |C 6 > 

We now compute the TTN description of \Cq) w.r.t. 
the tree T depicted in Fig. ^ using the procedure out- 



lined in theorem U] Consider the following Schmidt de- 
compositions of \Ce): 



|C 6 > = 4)El^ 1) )^ld 1) >3456 (Al) 
= iEl^VlefV (A2) 



iEl^WI^W (A3) 



X 



These decompositions are taken w.r.t. the bi- 
partitions ({1,2}, {3, 4, 5, 6}), ({1,2, 3}, {4, 5, 6}) and 
({1, 2, 3, 4}, {5, 6}), respectively; these correspond to the 
bipartitions (Aj,,Bj,), where e runs over all inner edges 
of T. All Schmidt vectors in IjAlfl are normalized, and 
the x are the square roots of the Schmidt ranks of the 
corresponding bipartitions |3Cj . 

We now show how the TTN description of \Cq) w.r.t 
the tree T is obtained, by applying the procedure pre- 
sented in theorem First, note that the depth A of T is 
equal to 3. We start by considering the single inner ver- 
tex of depth 3; this is the vertex which has leaves 1 and 2 
as lower vertices. We then compute the Schmidt decom- 
position ijAlf) . corresponding to the bipartition which is 
obtained by deleting the upper edge of this vertex. In 
a second step, we consider the single vertex in T having 
depth 2, and compute the corresponding Schmidt decom- 
position (|A2|) . Moreover, we write 



\C 6 ) 



x (i) x (2; 



El^K^Vfmf}, (A4) 



(where we have omitted the subscripts of the Schmidt 
vectors). Finally, we consider the Schmidt decomposition 
(|A3(I (corresponding to the uper edge of the unique depth 
1 vertex), and write it as 



A u 



(A5) 



Combining eqs. (|A5|) and l|A4(l then shows that \C&) can 
be written as follows 1311: 



1 



A A A ijk 

(A6) 

Note that the states (^j 1 ^ \4>f^) are defined on qubit 3, for 

(3) (2) 

every i and j, and that the states \Q ) are defined 
on qubit 4, for every j and k. 

Next we explicitly compute the Schmidt coefficients 
and Schmidt vectors in the above expansions, using the 
stabilizer formalism. 
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As for the Schmidt coefficients, note that 

2 = rank F2 r({l, 2}, {3, 4, 5, 6}) 
= rank F2 r({l,2,3},{4,5,6}) 
= rank F2 r({l,2,3,4},{5,6}), (A7) 

and therefore (using J7J)) all the Schmidt ranks of the 
above bipartitions are equal to 2 2 = 4. Thus, the indices 
k in eq. (|A6I) all run from 1 to 4, and we also have 



X (D =X (2) =X (3) = V4 = 2. 



(A8) 



It will be convenient to write the indices i, j, k as pairs of 
bits, and we will use the notations i = ab, j = cd, k = e/, 
where a, b,c,d £ {0, 1}. 

We now consider the Schmidt vectors w.r.t. the 
above bipartitions. We start with the bipartition 
({1, 2}, {3, 4, 5, 6}). Here one finds that 



Tr {3 ,4,5,6}(|C 6 )<C 6 |) = \l. 



(A9) 



Thus, a Schmidt basis for the subset {1,2} could simply 
be chosen to be the computational basis; in other words, 
we take 



\<j>§) = \a)®\b) = \ab), 



(A10) 



defined on the qubits {1, 2}, for every a, b £ {0, 1} . 

The same argument can be repeated for the vectors 
{1^'}, where we can take l^e/*) = l e /)> defined on the 
qubits {5, 6}, for every e,f£ {0, 1}. 

As for the bipartition ({1, 2, 3}, {4, 5, 6}), one can eas- 
ily show that 

Tr {4 ,5,6}(|C 6 )(C* 6 |) = -(/ + a z ® <r x <g> a z ) (All) 

and that, hence, the states 

\4®) = °Z®I®oi\Ls) (A12) 

form a valid Schmidt basis, where c,d £ {0, 1} and where 
I-L3) is the linear cluster state on 3 qubits, defined on the 
qubits {1,2,3}. 

To compute the vectors {\Q d )}, note that one has 



A( 2 ) 



cd \ C e 



(A13) 



Therefore, we have to compute expressions of the form 

[«L 3 |^OJ(8)^)®J]|C 6 ), (A14) 

for every c, d — 0, 1. To do so, we use that every n-qubit 
graph state |G) with adjacency T can be written as Q 



\G) = ^- 2 E (- 1 ) 



9g(«) I 



(A15) 



where {\u) \ u £ {0,1}"} is the n-qubit computational 
basis and where 



(A16) 



One then finds that (| A 14(1 is equal to (omitting multi- 
plicative constants) 

j ^_^qc 6 (x,y,z,u : v,w)+q L3 {x,y,z)+xc+zd j \ uvw y 

u,v,w \x,y,z / 

(A17) 

Straightforward algebra then shows that the power of — 1 
in the above expression is equal to 

x(w + c) + z(d + u) + qL 3 (u, v, w). (A18) 

Moreover, one has 



2 3 w — c and d = 
else 



^ ^ (_ ^x(w+c)+z(d+u) _ _f 



We then find that (|A17|I is equal to 

id) m-ir^V) j \c), 



(A19) 



(A20) 



for every c, d = Q, 1. Thus, these 4 states form the set 

{|f2?)}*=l, denned on qubits {4, 5,6}. 

The only remaining task is the computation of the 

states (4>^ \4>^cd) an( i (£e/ IC?)- To compute the former 
of these states, it follows from the above that one has to 
compute, for every a, 6, c, d £ {0, 1}, overlaps of the form 

W'iMrf) = ((a\®(b\®I)(a c z <S>I®a d z \L 3 )) 

= (-l) QC (a| ® (b\ ® o*\L 3 ). (A21) 

Using the expansion IjAlSjl . it is then easy to show that 
(IA21|I is equal to 



(~l) ac 'y^^—l) qL 3 ( - a ' b ' v ^ +dv \v) : 



(A22) 



for every a,b,c,d £ {0,1}, and these states are defined 
on qubit 3. A similar calculation can be performed to 
obtain 



(&\^) = 8fA-l) 9Lsid ' e ' c) \d), 



(A23) 



for every c,d,e,f £ {0, 1}, and these states are defined 
on qubit 4. 

We can now write down the TTN description of the 
state \Cq) w.r.t. the tree T depicted in Fig. ^ 



c *) = ^ E {1^)12 (ec- 1 

abcdef K \ v 



^ac+q L3 (a,b,v)+dv^ v ^ j x 



«G{0,1}" 



(Sf, c (-iy^ d ^\dU) |e/) 56 }, 

(A24) 
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where we have again indicated subscripts to specify on 
which qubits the states are defined. 
Recalling the definition of q^ 3 , namely 

<lL 3 {ti,t 2 ,h) ■= hh +t 2 h, (A25) 

for every ^,£2,^3 S {0,1}, we recover expression (|59ll . 
Note that one can easily check that Ij59(l is correct, by 
summing out all indices a, . . . , /: 

( Xl ...x 6 \c 6 ) = 1 J2 {K^ X2 (-^) ac+ab+bX3+dX3 

abcdef 

_(_-! \a;exi+xiX2-\ \-xs%6 

= _(-l)9C 6 (xi,..., :C 6) ) ( A26 ) 



where in the last equality we indeed obtain the correct 
computational basis expansion of |Cg). 
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